Boundary and impurity effects on entanglement of Heisenberg chains 
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We study entanglement of a pair of qubits and the bipartite entanglement between the pair and 
the rest within open-ended Heisenberg XXX and XY models. The open boundary condition leads 
to strong oscillations of entanglements with a two-site period, and the two kinds of entanglements 
are 180 degree out of phase with each other. The mean pairwise entanglement and ground-state 
energy per site in the XXX model are found to be proportional to each other. We study the 
effects of a single bulk impurity on entanglement, and find that there exists threshold values of the 
relative coupling strength between the impurity and its nearest neighbours, after which the impurity 
becomes pairwise entangled with its nearest neighbours. 
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I. INTRODUCTION 

Over the past few years much effort has been put into 
studying quantum entanglement in various quantum spin 
models. At zero temperature, entanglement naturally ex- 
ists in many-body ground states. The systems contain 
two classesjinite spin clusters [J0SSI1I1HEHII1 
ElEEiQHEOEl such as a ring of qubits, and 
infinite spin systems [m l2flj^. El 151 Im El l2r| where 
quantum phase transition |27j may occur. For a ring 
of qubits interacting via the Heisenberg XXX Hamilto- 
nian, the concurrence for the ground state, quantifying 
entanglement of a pair of qubits, is equal to the absolute 
value of the ground-state energy per site [1 . Universal 
scaling behaviours of quantum entanglement of ground 
states emerge at a quantum phase transition point in the 
anisotropic XY model [Tl I2I l2l| . Experimentally, en- 
tangled state of magnetic dipoles has been found be to 
crucial to describing magnetic behaviours in a quantum 
spin system pgj . 

In most of the previous studies on entanglement of 
many-body ground states, a periodic boundary condition 
(PBC) is assumed for spin chains. In compounds such as 
Mni2 and Fes, the metal ions within a single molecule 
form almost a perfect ring, which can be described by 
the Heisenberg interaction [2!j. On the other hand, spin 
chains with an open boundary condition (OBC) have 
been used to construct spin cluster qubits [H 0H| f° r 
quantum computation and employed for quantum com- 
munication from one end to another |32l l3^ |. Perfect 
state transfer have been obtained via the open chain 
without requiring qubit coupling to be switched on and 
off |33j • These investigations reveal that open chains are 
of great advantage in implementing quantum information 
tasks. Thus, the study of the entanglement structure 
in open spin chains will be of importance as the entan- 
glement underlies operations of quantum computing and 
quantum information processing. Here, we investigate 
the open boundary effects on ground-state entanglement. 

A linear open A-qubit chain can be viewed as a ring 
of (A + l)-qubit chains with an impurity ((A + l)-th 



qubit), where the coupling <//v+i,i = Jn+i.n = 0. The 
impurity plays a very important role in condensed matter 
physics, and it will be interesting see its effect on quan- 
tum entanglement. In earlier studies, impurity effects on 
entanglement have been considered with a small three- 
spin system |34j . Here, we consider large spin systems 
described by the Heisenberg XY model. 

The paper is organized as follows. In Sec. IIA, by exact 
diagonlization method, we study the ground-state entan- 
glement in the Heisenberg XXX model with an OBC. 
The mean nearest-neighbour entanglement is found be 
proportional to the ground-state energy per site. For 
studying systems with large number of site, in Sec. IIB, 
we consider the XY model with an OBC, which can be 
solved exactly by Jordan- Wigner mapping [35l |. Entan- 
glements display strong oscillations over lattice. With 
increase of the ground-state energy, the mean pairwise 
entanglement increases in the XY model, while decreases 
in the XXX model. In Sec. Ill, we study impurity ef- 
fects on entanglement of the spin systems described by 
the XY model, where the impurity is located in the bulk 
of the system. We conclude in Sec. IV. 



II. BOUNDARY EFFECTS 



A. Heisenberg XXX model 

We consider a physical model of A qubits interacting 
via the isotropic Heisenberg Hamiltonian with an OBC, 

HxXX — — ^ (l+Crix&i+lx+&iyO'i+ly+0'izC r i+lz) , (1) 
i=l 

where (7j = (cr^ , a iy , a iz ) is the vector of Pauli matrices 
and J is the exchange constants. The positive and neg- 
ative J correspond to the antiferromagnetic (AFM) and 
ferromagnetic (FM) case, respectively. 

The SU(2) symmetry ([H,S a ] = 0) is evident, where 
S a = a ia/%, a £ { x , V, This symmetry guaran- 

tees that reduced density matrix py of two qubits, say 



2 



qubit i and j, for the ground state p has the form 



Pij 



I mj 



Vo 







(2) 







in the standard basis {|00), |01), |10), 1 11)}. 

From the reduced density matrix, the concurrence_|3£ 
quantifying the entanglement is readily obtained as 



Cij =2 max [0, \z^ \ - 

= max[0,|G?j|-G?5/2-l/2] 
= max(0,|Gf/|-Gf//2-l/2) 
= max(0,-3G??/2-l/2), 



(3) 



where Gff 



Tic(ai a (J2ap) are correlation functions. The 
second equality follows from the relations between matrix 
elements of and correlation functions given by = 
(1 + Gf/)/4 and z tj = Gff/2. The SU(2) symmetry lead 
to the third equality. The fourth equality follows from 
the inequality \Gfj \ < 1, which is a special case of a more 
general result that | (A) \ < 1 for any Hermitian operator 
A satisfying A 2 = 1. We see that the qubit i and j are 
entangled if G?J < -1/3. 

For studying bipartite entanglement, we consider even 
number of sites, for which ground states of the Heisenberg 
model are nondegenerate, and thus pure. The bipartite 
entanglement of pure states is well-defined by entropies 
of one subsystem. Here, we choose the linear entropy 
to quantify the bipartite entanglement between a pair of 
qubits and the rest. From Eq. |(SJ|, in a similar way to 
get Eq. J3J), the linear entropy for /ay is obtained as 



Ej 



1 



Tr(p?) 



1 



l--[l + 3(Gf/) 2 ]. (4) 



The reduced density matrix pij is only determined by the 
correlation function Gff, so do the concurrence and the 
linear entropy. 

Due to the nearest-neighbour nature of the interac- 
tion in our system, the entanglement between a pair of 
nearest qubits is expected to be prominent comparing 
with a pair of non-nearest-neighbour qubits. Thus, we 
focus on the nearest-neighbour case in the following dis- 
cussions. For the isotropic Heisenberg Hamitonian with 
a PBC, entanglements between qubit i and i + 1 are in- 
dependent on index i. However, for the case of OBC, 
the concurrence G^ + i and the linear entropy E ii+ i must 
be site-dependent due to the breaking of translational 
symmetry. Next, we study entanglements using exact 
diagonalization method. 

We first diagonalize Hamitonian Hxxx to obtain the 
ground state, from which we calculate the correlation 
function GS_ X , and then the pairwise and bipartite en- 
tanglements via Eqs. © and Q}, respectively. The nu- 
merical results for 12 qubits (4069-dimensional Hibert 
space) are shown in Fig. 1. The entanglements oscillate 




FIG. 1: Ground-state nearest-neighbour concurrence (square 
line) and linear entropy (star line) versus site number for L = 
12 in the Heisenberg XXX model. 



with a two-site period, and the concurrence and the lin- 
ear entropy are 180 degree out of phase with each other. 
The pair with qubits 1 and 2, and the one with qubits 
L— 1 and L display maximal pairwise entanglement (min- 
imal bipartite entanglement). For two qubits interacting 
via the Heisenberg Hamiltonian, they prefer to the max- 
imally entangled singlet state |^) = 1/V2(|01) - 1 10)) 
with concurrence C = 1. Qubit 1 favor maximal en- 
tanglement with the only nearest-neighbour qubit 2, and 
result in higher concurrence C\2- For qubit 2, there are a 
competition between qubit 1 and 3, and they both favor 
maximally entangled with qubit 2. Qubit 2 shares large 
entanglement with qubit 1, and thus results in less entan- 
glement with qubit 3, and the oscillatory feature appears. 
As the pair of qubits 1 and 2 have strong entanglement, 
the bipartite entanglement between the pair of qubits 
with the rest of the system is relatively suppressed as we 
have seen from the figure. 

It will be instructive to give two extreme examples 
to reveal relations between pairwise entanglement and 
the biparite entangelement between the pair and the 
rest. The first one is the singlet state with C — 1 and 
E = 0, and second one is the maximally mixed state 
p u = diag(l/4, 1/4, 1/4, 1/4) with G = and E = 0.75 
(maximal linear entropy). Thus, in general, the above 
numerical results and the two simple examples suggest 
that the more the pairwise entanglement, the less the 
bipartite entanglement. Of course, this is not a general 
statement for arbitrary multi-qubit states. 

Next, we consider the mean entanglement and study 
its relations with ground-state energy £q. The nearest- 
neighbour mean concurrence and linear entropy can be 
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FIG. 2: Mean concurrence (square line) and mean linear en- 
tropy (star line) versus L and eo in the Heisenberg XXX 
model. 



defined as 



L-l 



a 



E, 



t=i 

1 



(5) 



We made numerical calculations for even L from 4 to 12, 
and the entanglements versus L or eo are shown in Fig. 2, 
where eo = £b/£ is the ground-state energy per site. The 
pairwise entanglement decreases as L or eo increases, and 
so does the bipartite entanglement. More interestingly, 
the pairwise entanglement decreases linearly with eo- We 
may obtain an analytical relation between C m and eo as 
follows. 

The numerical results for the finite lattices show that 
Cm+1 > 0, and thus from Eq. we obtain 



(~1ZZ 



-2C; 



1 



(6) 



From Hxxx Q an d after taking into account the SU(2) 
symmetry, we get 



L 1 3 v -v z ^ 

i=i 



Substituting Eq. © into leads to 



eo 



L-l 
L 



C m (L), 



(7) 



(8) 



where we have used the definition of mean concurrence. 
The above relation shows that the mean entanglement is 
proportional to the ground-state energy per site, imply- 
ing that less energy gives more pairwise entanglement. 
Although this relation is obtained for finite lattices, we 



make a conjecture that it is valid for any number of lat- 
tice sites. 

The concurrence and linear entropy obtained for finite 
sites of lattice up to 12 are not so close to those for infinite 
lattice given by 



C m (oo) 
E m (oo) 



=2 In 2 
2 r . 



- 1 w 0.3863, 

In 2(2 In 2 - 1)] « 0.4882. 



(9) 



This is due to the finiteness of lattice, and also due to 
the effects of boundaries. To investigate effects of open 
boundaries on entanglement in larger systems, we next 
study quantum Heisenberg XY model, which can be 
solved exactly by the Jordan- Wigner transformation [35| . 



B. Heisenberg XY model 



by 



The Heisenberg XY Hamiltonian with an OBC is given 



Hxy — — 



2 ^ 



(10) 



There are two symmetries in the Hamiltonian. One is a 
U(l) symmetry {[H,S Z ] = 0), and another is a Zi sym- 
metry ([H, £J = 0), where S x = ai x ® a 2x ® • ■ • <8> &Nx- 
They guarantee that reduced density matrix py of two 
qubits i and j for the ground state p is given by Eq. J5J). 
For the lack of SU(2) symmetry, the correlation func- 
tions G** and G?? are no longer equal, and then the 
concurrence is given by 



C tf = max [0,10^1-0^/2 -1/2] 



(11) 



determined by two correlation functions. 

For nearest-neighbor cases, two correlation functions 
and Gff +1 are dependent, and the latter can be 
written in terms of the former as we will see shortly, 
and thus the concurrence is only determined by a single 
correction function Gff +1 . To show this, we use the well- 
know Jordan- Wigner mapping |35j . under which Hamil- 
tonian Hxy exactly maps to 



L-l 



Hxy =J^2 ( c I Cl +! + c l+i c i) 



(12) 



due to the OBC, where Ci and cj are fermionic annihi- 
lation and creation operators, respectively. Then, in the 
fermionic representation, the two correlation functions 
can be written as 



J, 



Zli+l ={cic i+ i), 

u ii+ i =(qtc i c| +1 c j +i) = 1/4- | (4c l+1 ) | 2 , (13) 

where we have used Wick theorem 37], and the fact 
(c\ci) = 1/2 for any i. Thus, the concurrence is de- 
termined only by one correlation function (c|ci+i), 



a 



ii+l 



2max(0, |(qcf 



1/4) (W) 



4 



To compute the correlation function (cjcj-|_i), we first 
diagonalize Hamitonian Hxy by the following transfor- 
mation 



k=l 



sin | - I c fc (15) 



k=l 



L + l 



After the transformation, Hamiltonian Hxy can be writ- 
ten as 

H *Y = J2<^kclc k = ^-2 cos fy^y] c{c k , (16) 



k=l 



fc=l 



where is the single-particle energy, and J = — 1 is 
assumed for convenience. The ground state is then ex- 
pressed as 



\*) GS = c\cl...,cl /2 \0), 



(17) 



with energy Eq — 5Jfc=i e k- From the above expression 
of the ground state, the correlation function is obtained 

as 



(4 c i+i> 



L/2 

/]9i,k9i+l,k- 
k=l 



(18) 



The combination of Eqs. 1|14|) and (|18|l gives exact ana- 
lytical expression for the concurrence. By using Eq. (|13|l , 
the linear entropy of the two-qubit reduced density ma- 
trix /Ojj+i is obtained as 



E u+1 = 3/4- 4|(c J t c, +1 )| 4 - 2\(4c l+1 )\ 2 



(19) 



which is also analytical. From these analytical expres- 
sions, entanglements are readily calculated numerically 
for a very large L. 

Figure 3 displays features of entanglement distribution 
over lattice sites, which are similar to those of the Heisen- 
berg XXX model. Strong oscillations of entanglements 
are more evident. The similarity arises since two qubits 
interacting via the Heisenberg XY interaction also favor 
the singlet state. Near the bulk area, the entanglements 
oscillate with nearly the same amplitude with respect to 
the mean value, and the boundary effects diminish. 

Now, we study the mean entanglements for different L. 
Figure 4 gives the entanglement versus L and the energy 
per site. With the increase of L, the mean concurrence 
decreases, while the linear entropy increases, and both 
saturate for large L. For a infinite lattice, the concur- 
rence and the linear entropy are given by 



C ro (oo) =- + \ - ~ pb 0.3393, 

E m (oo) =^ - 4 - 4 w 0.5063. 
4 7r 4 n z 



(20) 



We see that the concurrence and linear entropy with L = 
200 are nearly identical to those with infinite L. From 
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FIG. 3: Ground-state nearest-neighbour concurrence (square 
line) and linear entropy (circle line) entanglement versus site 
number for L — 100 in the Heisenberg XY model. 




FIG. 4: Mean concurrence (square line) and mean linear en- 
tropy (star line) versus L and eo in the Heisenberg XY model. 



Fig. 4, we observe that the mean concurrence increases, 
while the mean linear entropy decreases with the increase 
of energy. The behaviour of the mean concurrence is 
opposite to that in the Heisenberg XXX model. Here, 
the lower energy does not correspond to higher pairwise 
entanglement. 



III. IMPURITY EFFECTS 

As we have stated in the introduction, a linear open 
iV-qubit chain can be viewed as a ring of (N + l)-qubit 
chains with an impurity. In this section, we study single 
impurity effects on entanglement of Heisenberg chains. 
We consider the Heisenberg XY model described by the 
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Hamiltonian 

*f-i j 

HxY = ^ — ^ (cix&i+lx + Piy&i+ly) ■ (21) 

i=l 

with an OBC. Now, we assume the single impurity spin 
is located in the bulk (site N/2) j3^|, namely, 

Jn/2-1,N/2 — Jn/2, N/2+1 = J = a J, 

= J = 1 for other i, (22) 

where a characterizes the relative strength of the cou- 
pling between the impurity qubit and its nearest neigh- 
bours. To be specific, we focus on pairwise entanglement 
in the following discussions. 

Due to the existence of a single impurity in the bulk, 
we do not expect analytical results for entanglement. 
However, the Hamiltonian still have the U(l) and Z 2 
symmetries, we can map the Hamiltonian exactly to a 
fermionic one, and the Wick theorem applys. Thus, 
Eqs. (|llfl and (|14fl for the concurrence are valid and ap- 
plicable to the impurity model. All we need to is to 
diagonalize a L x L matrix to obtain the coefficients 
g n k in the expression c n — X}fc=i SnfcCfc) and then cal- 
culate the correlation function (cjcj+i) via the relation 

(c|ci_|_i) = J2k=i 9ik9i+i,k- Note that the summation is 
from k = 1 to L/2 since the system considered is non- 
degenerate (the single-particle energy e& is not zero) and 
number of the negative values of e& are L/2. 

From Fig. 5, we see that the impurity leads to addi- 
tional oscillations of the pariwise entanglement in the 
bulk region. For small a = 0.1, the concurrences 
C l A r /2-i.JV/2 an d Cat/2, N/2+1 are zero. Even the coupling 
between impurity qubit and its nearest neighbours are 
not zero, the entanglement between them vanishes due 
to competition among qubits. For larger a = 0.6, the en- 
tanglements between impurity qubit and qubits L/2 + 1 
and L/2 - 1 build up, where C N / 2 _ ltN / 2 > Cn/2,n/2-u 
which also holds for a = 1.2 > 1 and a = 2 > 1. The 
difference between the two concurrences results from the 
choice of even TV, leading to a non-symmetry of the en- 
tanglement distribution. We see that even one single 
impurity have strong effects on entanglement structure, 
especially in the region near the impurity. 

From the above analyses, it is expected there exists 
a threshold value of cv, after which the impurity qubit 
and its nearest neighbors becomes entangled. In Fig. 6, 
we plot the nearest-neighbour concurrences in the bulk 
region as a function of a. It is evident that there exist 
threshold values ol\ for concurrence /2-\.n /2 an( i a 2 
for concurrence C^/ 2 .n/2+\- The threshold value a 2 is 
slightly larger than ot\, and the pairwise entanglement 
between the impurity and qubit N /2 — 1 is always a lit- 
tle stronger than that between the impurity and qubit 
N 1 2 + 1 when a > a 2 . We also observe that the increase 
of a suppresses the concurrence Cm/2+i, N/2+21 while en- 
hance the concurrence Cm/2+2. n/2+z- When a = 1 (no 
impurity case), all concurrences shown in the figure are 
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FIG. 5: Ground-state nearest-neighbour concurrence versus 
site number for L = 100 in the Heisenberg XY model with 
an impurity. 




a 



FIG. 6: Ground-state nearest-neighbour concurrences 
Cjv/2-i,a72 (square line), C N/2 , n/2+i (circle line), 
C]v/2+i,jv/2+2 (diamond line), and C n/2 + 2 .n/2+3 (star 
line) as a function of a for the impurity model. 

almost identical since we have chose a large L = 500, 
which diminishes the boundary effects in the bulk and 
the amplitudes of oscillations become very small. 



IV. CONCLUSION 

In conclusion, we have studied ground-state entangle- 
ments in Heisenberg XXX and XY models with an 
OBC. The OBC leads to strong oscillations with a two- 
site period of entanglement in a pair of nearest-neighbour 
qubits and bipartite entanglement between the pair and 
the rest. The maximal pairwise entanglement and min- 
imal bipartite entanglement occurs at open ends, and 
the two kinds of entanglements are 180 degree out of 
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phase with each other. In both models, the two-qubit 
reduced density matrix is determined by only one corre- 
lation function, and so do the entanglements. We have 
found that the mean entanglement is proportional to the 
ground-state energy per site in the XXX model. With 
increase of the ground-state energy, the mean pairwise 
entanglement decreases in the XXX model, while in- 
creases in the XY model. 

We study the effects of a single bulk impurity on en- 
tanglement, and find that the impurity leads to addi- 
tional oscillations of entanglement in the bulk region. We 
also find that there exists threshold values of the relative 
coupling strength between the impurity and its nearest 
neighbours, after which the impurity becomes entangled 
with its nearest neighbours. As the entanglement under- 



lies operations of quantum computation and quantum 
information processing, the structures of entanglement 
found in the present studies are useful when we make 
a simulation of quantum systems where boundary and 
impurity effects cannot be negligible. 
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